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ABSTRACT 

We present a model for the CO molecular line emission from high redshift Submillimeter 
Galaxies (SMGs). By combining hydrodynamic simulations of gas rich galaxy mergers with 
the polychromatic radiative transfer code, SUNRISE, and the 3D non-LTE molecular line ra- 
diative transfer code, TURTLEBEACH, we show that if SMGs are typically a transient phase of 
major mergers, their observed compact CO spatial extents, broad line widths, and high excita- 
tion conditions (CO SED) are naturally explained. In this sense, SMGs can be understood as 
scaled-up analogs to local ULIRGs. We utilize these models to investigate the usage of CO as 
an indicator of physical conditions. We find that care must be taken when applying standard 
techniques. The usage of CO line widths as a dynamical mass estimator from SMGs can pos- 
sibly overestimate the true enclosed mass by a factor ~1 .5-2. At the same time, assumptions 
of line ratios of unity from CO J=3-2 (and higher lying lines) to CO (J=l-0) will oftentimes 
lead to underestimates of the inferred gas mass. We provide tests for these models by outlin- 
ing predictions for experiments which are imminently feasible with the current generation of 
bolometer arrays and radio-wave spectrometers. 

Key words: cosmology :theory-galaxies:formation-galaxies:high-redshift- 

galaxies:starburst-galaxies:ISM-galaxies:ISM-ISM:molecules 



1 INTRODUCTION 

Understanding the origin and evolution of active galaxy popu- 
lations at cosmological redshifts remains an outstanding prob- 
lem. Of particular interest are a population of submillimeter- 
luminous galaxies discovered via deep, blind s urveys with SCUBA 
on the JCMT at a median redshift of z ~ 2 dBarger et all 1 19981 : 
iHughes et all 1 19981 ; Ichapman et afl l2003al) . These Submillime- 
ter Galaxies (SMGs) are emp irically defined w ith the flux limit 
> 5-6 mj£] (e.g. iBlain et alj l2002h . and appear to 
be a highly clustered population of galaxies forming s t ars at 
prodi g ious rates (SFR > 1 3 M^yr 



Blain et al.1 



20041: ISwinbank et alj 12004 tM enendez- Deimestre et al] 



Kovacs et al 



2002, 



2007; 



20061 ; IValiante et alj|2007l ; ICoppin et alj|2008ah . In 



addition to undergoing prodigious star format ion, some SMGs 
additionally host heav ily obscured AGN (e.g. Iyison et all 1 20021 ; 
I Alexander et al]|2005al lb[ 120081 ; iBorvs et alj|2005h . These galaxies 
may contribute a substantial fraction of the cosmic star formation 
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density telainetalj[l99^.l2002h . as well as serve as prime candi- 
dates for studying the coevolution of black hole growth and star 
formation at an epoch of heightened galaxy formation and evolu- 
tion. 

Many open questions remain regarding the physical nature of 
SMGs. Are they isolated galaxies, or galaxy mergers? What is their 
potential place in an evolutionary sequence? Are they high-redshift 
analogs to local infrared luminous galaxies (e.g. ULIRGs)? What 
are their typical gaseous and stellar masses? 

Emission from the star forming molecular interstellar medium 
(ISM) has the potential to elucidate some of these questions. For 
example, 12 CO (J=l-0; hereafter, CO) emission can serve as a 
measure both for the total amount of H2 molecular gas avail- 
able for future star formation (e.g. iDownes & Solomonl 1 19981 ; 
Solomon & Barret tl ll99ll; iGreve et al.l 12005), as well as previous 
stellar mass assembly by serving as a dynamical mass tracer (e.g. 



stellar mass assembly by serving as a dynamical mass tracer (e.g. 
Tacconi et all2006l : iBouche et alll200l [lloll2007l ; iNaravanan et all 



2008b). The high resolution images available from millimeter- 
wave interferometry can reveal source sizes and morphologies (e.g . 
Genzel et alj|2Qol;lTaccom et al.ll2006l . 120081; I Younger et al]|2008l; 



Iono et al l l2009h . Additionally, the rotational ladder from H2 trac- 
ers such as CO reveal the ther mal conditions and me an densi- 
ties of the star forming ISM iWeifi et alj [2001 l2007t) . CO has 
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been d etected routinely in high r e dshift SMGs for o ver a decade 



(e.g. iFraver et all 1 19981 1 19991; iNeri et al 



2003; iGreve et al.ll2005l : iTacconi et alj|200d 



.20031 ; iGenzel et al.1 
2008, and references 



therein). During this time, a number of pioneering papers have af- 
forded the community a wealth of CO data from z ~ 2 SMGs. The 
general molecular emission properties from SMGs can be summa- 
rized as follows. 

Hi Gas Masses: SMGs are incredibly gas rich, with massive 
H2 gas reservoirs. Extrapolation from typical CO (J=3-2) rest frame 
measurements from SMGs suggests H2 masses of order 10 10 -10 n 
Mq, making these s ome of the most molecular gas rich galax- 
ies in the Universe (Greve et al. 2005; Solomon & Vanden Bout 
2004 ITacconi et alj I2OO6I ; ICarilli & Wand 120061 ; ITacconi etaU 



2008; Coppin et a l. 2008b). An analysis of semi-analytic models of 
SMG formation by Swinbank et al.l J2OO8I) found good agreement 



between the cold gas mass in their simulations and the observed 
values. 

Line Widths: The CO line widths in SMGs are typically broad 
with median line w i dth ran ging from ~600 - 800 km s" 1 (F WHM; 
Greve et alj |2005t ICarilli & Wand 120061; ITacconi et alj 120061 ; 



Coppi n et alj|2008bl ; ITacconi et alj|2008l ; llono et alj|2009l) . Com- 



parisons to CO-detected quasars at comparable re dshifts have had 
conflic ting results. Analysis of literature data by ICarilli & W ang 
d2006h has suggested that SMGs typically have broader line widths 
than z ~ 2 quasars by a factor of ~2.5 with KS tests indicat- 
ing the two are not drawn from the same parent population at 
the >99% confidence l evel. On the other hand, new detections 
bvlCo ppin etaD |2008b) have found rather similar CO line widths 
from quasars and SMGs, and KS tests showing that the two arise 
from the same parent population at the >95% confidence level. 

Images: Interferometric CO imaging of SMGs has shown 
the spatial distributio n of most SMGs to be relatively compact 
(Rhwhm < 1-2 kpc; ITacconi et al.ll2006l) . Some individual so urces 
have shown rather extended emission (e.g. lGenzel et al .120031) . Cu- 
riously, images have shown ev idence for compact disk-like mo- 
tion i n at least a few sources jGenzel et al.1 120031 ; ITacconi et al.l 
2006), as well as extende d emission in what appears to be interac t- 
ing/merging galaxies (e.g. Ilvison et alJl200ll ; ITacconi et alJl2006h . 

Excitation Properties: While very few SMGs have been ob- 
served in multiple CO lines, the existing data suggests that these 
galaxies exhibit relatively high molecular excitation conditions. 
The CO line spectral energy distributions (CO SEDs; alternatively 
known as CO rotational ladders) are see n to typically turn over 
at the J=5 level ( ~83 K above groun d; IWeifi et alj|2005l , 120071 ; 
iGreve et alj2005l ; lHainline et al.l200"6l) . That said, fe w objects have 
been detected in their rest-frame CO (J=l-0) line Mainline et al.1 
2006). As H2 masses are typically inferred from CO (J=l-0) mea- 
surements, the excitation properties of higher lying levels (e.g. line 
ratios) with respect to the ground state transition are crucial to con- 
strain for the purposes of deriving molecular gas masses. 

While the information provided by CO observations of 2 ~ 2 
SMGs is indeed invaluable, questions remain regarding the in- 
terpretation of many of the various aforementioned observational 
characteristics of these sources. For example, what is the true rela- 
tionship between CO line widths from SMGs and quasars? Is there 
an evolution in the CO line widths of SMGs? How reliably can the 
CO line width be used as dynamical mass tracer in these galaxies? 
How can observations of higher lying CO lines (in the rest frame) 
be extrapolated to the rest fra me CO (J=l-0) lumin osity in order 
to derive an H2 gas mass (e.g. Hain line et ai] |2006)? Can the ob- 
served CO emission line widths, excitation properties, and images 
be explained by a merger-driven scenario for SMG formation and 



evolution (e.g. "scaled up ULIRGs"; lTacconi et al.ll2006ll2008h ? It 
is clear that a theoretical interpretation behind the molecular line 
emission properties may be valuable for elucidating some of the 
aforementioned issues, as well as providing interpretation for forth- 
coming observations. 

Along with providing interpretation for observations, theoret- 
ical calculations of CO emission from simulated SMGs can pro- 
vide direct tests of mo dels of SMG formation and evolution. In 
iNaravanan et alj d2009h . we presented a merger-driven model for 
the formation of SMGs which reproduced the full range of ob- 
served 850 \im fluxes from SMGs, the optical-mm wave SED, and 
characteristic stellar, black hole, and dark matter masses. Compar- 
ing the simulated molecular gas properties of these model SMGs 
to the extensive data sets in the literature provides a strict test of 
the models. In this paper, we investigate the CO emission proper- 
ties from SMGs by combining these SMG formation models with 
3D non-local thermodynami c equilibrium (LTE) molecular li ne ra- 
diative transfer calculations l lNaravanan et a l. 2006a, 2008eJ). The 
goals of this paper are to: (a) provide direct tests of merger-driven 
formation mechanisms for S MGs by comparin g the s imulated CO 
emission from the models of INaravanan et all |2009) to observa- 
tions, and, given a sufficient correspondence between models and 
observations, (b) provide interpretation for existing and future ob- 
servational data. 

This paper is organized as follows: in § [2] we describe our 
numerical methods for our hydrodynamic, molecular line radiative 
transfer, and polychromatic SED radiative transfer simulations. In 
§[3] we summarize the evolution of the submillimeter, 5-band, and 
H2 properties of our simulated galaxies. In W. II we study the CO 
morphology, molecular disk formation and CO spatial extents in 
SMGs. In § 14.21 we explore the origin of the broad observed CO 
line widths. We use these results to investigate the usage of CO as 
a dynamical mass tracer in § [5] and analyze the CO excitation in 
SMGs in § [6] In § [7] we provide imminently testable observational 
predictions, and in §[8] we summarize. 



2 NUMERICAL METHODS 

Generically, our methodology involves three steps. First, we simu- 
late the hydrodynamic evolution of galaxies utilizi ng the smoothed - 
particle hydrodynamics (SPH) code, GADGET- 3 JSpringelll2005l) . 
We then calculate the molecular line emission properties and SEDs 
from these simulated galaxies in post -processing utilizing th e ra- 
diative tra nsfer codes TURTLE BEACH (Naravanan et al. 20083) an d 
SUNRISE Jjonsson et alj|2009l) . This method is summarized in the 
flowchart presented in Figure [T] In this section, we present the de- 
tails of the hydrodynamic and radiative transfer methods. 



2.1 Hydrodynamics 

A major goal of our program is to understand the evolution of 
model SMGs while remaining constrained to observations for phys- 
ical input to the simulations. When data for SMGs is unavailable, 
we turn to observational constraints from the Galaxy or local star- 
bursts. We consider the formation of SMGs in gas rich binary 
galaxy mergers at high redshift. This is motivated by observed ra- 
dio, CO, and optical morph ologies of SMGs which a ppear to show 
signs of in teractions (e.g. IChapman etai] l2003bl : ITacconi et al.1 
2006, 20081) , as well as theoretical models which demonstrate that 
mergers serve as an efficient means of triggering nuclear starbursts 
tearnes & Hernquist|[l99ll 1 19961 : iMihos & Hernquistlll994 1 19961 : 
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Table 1. SMGs are ordered with decreasing halo mass. Column 1 is the name of the model 
used in this work. Columns 2 & 3 are initial orientations for disk 1, Columns 4 & 5 are for 
disk 2. Column 6 gives the virial velocity of the progenitors, Colu mn 7 gives their halo 
masses, and Column 8 the mass ratio of the merger. Please see iNaravanan et alj 120091) 
for generalized results regarding the physical properties of these galaxies. 



Model 


6>i 


<t>i 


02 


4>2 


V c 
(km s" 1 ) 


M 


Mass Ratio 


SMG1 


30 


60 


-30 


45 


500:500 


2.5xl0 13 


1:1 


SMG2 


360 


60 


150 





500:500 


2.5 xlO 13 


1:1 


SMG3 


-109 


-30 


71 


-30 


500:500 


2.5xl0 13 


1:1 


SMG4 


30 


60 


-30 


45 


500:320 


1.6xl0 13 


1:3 


SMG5 


360 


60 


150 





500:320 


1.6xl0 13 


1:3 


SMG6 


-109 


-30 


71 


-30 


500:320 


1.6xl0 13 


1:3 


SMG7 


30 


60 


-30 


45 


500:225 


1.4xl0 13 


1:12 


SMG8 


360 


60 


150 





500:225 


1.4xl0 13 


1:12 


SMG9 


-109 


-30 


71 


-30 


500:225 


1.4xl0 13 


1:12 


SMG10 


30 


60 


-30 


45 


320:320 


6.6xl0 12 


1:1 


SMG11 


360 


60 


150 





320:320 


6.6xl0 12 


1:1 


SMG12 


-109 


-30 


71 


-30 


320:320 


6.6xl0 12 


1:1 


SMG13 


30 


60 


-30 


45 


225:225 


2.35 xlO 12 


1:1 



GADGET 



SPH simulations of galaxy 
evolution to calculate density, 
temperature and velocity information 
of high redshift galaxies in evolution. 




TURTLEBEACH 

3D non-LTE molecular line radiative 
transfer. Returns molecular line 
spectra, maps and excitation 
properties of model galaxies simulated 
with GADGET. 



SUNRISE 

3D polychromatic dust radiative transfer. 
Calculates transfer of UV-mm wave 
radiation from stars and AGN through 
dusty ISM and returns spatially resolved 
SED information. Used here 
to determine when galaxies appear 
as SMGs. 



Figure 1. Flowchart summarizing methodology. The ga laxies are simu - 
lated hydrodynamically using the SPH code, GADGET-3 jSpringeil2005h . 
The simulated molecular line emission and SEDs are then calculated 
in post-processing utili zing the radiat ive transfer codes TURTLEBEACH 
I Naray anan et alj2008al) and SUNRISE fjonsson et alj20 09). respectively. 



Sprin gefet al.l l2005h . Moreover, simulations by INaravanan et all 
( 120091) have suggested that isolated galaxies and minor mergers 
(mass ratio < < 1 : 10) are unlikely to result in a Ssso >5 mjy SMG. 

The hydrodynamic simulations were performed with the 
smoo thed-particle hydrodynamics code, GADGET- £| dSpringell 
120051) which utilizes a fully conservative SPH formalism 
( Sprin gel & Hemquistl |2002|) . The hydrodynamic sim ulations in- 
clude prescriptions for radiative cooling of the gas dKatz et all 
Il996l : iDave et alj 1 19991) . and a multiphase ISM in which cold 



The main improvement in GADGET-3 over GADGET- 2 (described by 
ISpringelll2005l) is better load balancing on parallel processors. 



clouds are considered to be in pressure equilibrium with hot gas 
( McK ee & OstrikeJ 19771 : ISpringel & Hernquisfe oOB ) . This multi- 
phase ISM is implemented such that cold clouds grow through ra- 
diative cooling of hot gas, and heating fro m star formation can 
evaporate cold clouds 1 Sp ringel et al.ll2005r) . The effect of super- 
novae on the ISM is treated via an effective equation of state (EOS). 
Here, we employ the full multi-phase EOS where supernovae- 
driven pressure optimally maintain disk-stabilit y (gEOS=l; For 
more details, see Figure 4 of ISpringel et ai1l2005l) . 

Star formation occurs following a volum etric generalizati on 
of the Kennicutt-Schmidt relation, SFR oc n 1 5 dKennicuttll998dlbl : 
iKennicutt etai]|2007l : ISchrmddfl959l: ISpringel & Hernquistll2003h . 
which results in disks consiste nt with the local surface-density scal- 
ing relations dCox et alj2006h . The star formation timescale is cho- 
sen such that isolated disk models are consistent with the normal- 
ization of the local Kennicutt-Schmidt relations. While there are 
rather few measurements of observed star formation rate (SFR) 
relations at cosmological redshifts, tentative evidence exists that 
z ~ 2 galaxies may lie on the present-day Kennicutt-Schmidt rela- 
tion (Bou che et ai]|2007l) . Indeed, this is theoretically favorable as 
the relation SFR~n 1,5 may be understood in terms of free-fall time 
arguments which are redshift invariant. 

Black holes are included in the simulations as sink particles 
which accrete following a Bondi-Hoyle-Lyttleton parameterization 
according to the lo cal gas density and sound speed (Bondi & Hovle 
ll944l : lBondilll952l) . The bolometric luminosity of the black hole is 
set at L\, \=eM<? , where e=0. 1. Feedback from the black hole(s) 
is modeled such that a fixed fraction of this luminosity (here, 5%) 
couples thermally and isotropically to the surrounding ISM. This 
coupling efficie ncy is set to match the normalization of th e local 
M - a relation dDi Matteo et alj|2005l : lspringel et al.1l2005l) . 

The progenitor disk galaxies are initialized wi th alHernquist 

(1990) dark matter halo profile (for details, see Sprin gel et al. 
2005) with virial properties scaled to be appropriate for z ~ 3 (such 
that th ey may represent galaxies undergoing a major merger by 
2 ~ 2: lBullock et al.l200ll : lRobertson e t al. 2006). The progenitors 
begin with an initial gas fraction of / fl =0.8, resulting in a gas frac- 
tion of / g ~ 20-40% by the time the galaxies approac h final coales- 
cence. This is comparable to recent measurements bv lBouche et alj 
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d2007l) and lTacconi et al, I d2008l) who find tentative evidence for gas 

fractions in z ~ 2 SMGs of ~ 40%. 

Here, we consider the iNaravanan et alj 12009) series of 12 
models, varying total mass, mass ratio, and orbit. We include one 
additional lower mass merger which will not make an SMG (model 
SMG13), though will be useful for predictions of lower flux galax- 
ies (§ |7]l. The initial conditions of the merger simulations are 
summarized in Table [2~T1 The primary progenitors are initialized 
with circular velocities ranging f rom 1^=320-500 km s~ , simi- 
lar to measurements tabulated by iTacconi et alj d2008l) . This cor- 
responds to halo masses of order Mdm ~10 12 -10 13 M Qj _ com- 
parable to those inferred by clustering measurements (Bl ain et al.1 
12004 IS winbank et alj|2008l) . We model the formation of SMGs in 
1:1 mergers in ~10 12 , 5x 10 12 , and ~1O 13 M halos, and 1:3 and 
1:12 mergers with a Mdm ~ 10 13 Mq primary galaxy. We uti- 
lize 3 different initial orbits for the galaxies. For clarity, throughout 
this work, we primarily focus on the results from a fiducial merger 
(in Table 12. ll model SMG10) which produces a relatively average 
S 85 o«5-7 mJy SMC0. The CO results (next section) are generic 
for all models considered here. There is a dispersion amongst the 
simulated submillimeter fluxes when v arying merger orbit (eve n at 
a constant galaxy mass; see Figure 2 of lNaravanan et al.1 ( 120091 ) for 
a direct quantification of this dispersion). The fiducial SMGs stud- 
ied here lie in the middle of th is dispersion for g iven h alo masses, 
and, as such, are typical. See INaravanan et ak I d2009h for results 
regarding the physical properties of these merger simulations. 



2.2 Radiative Transfer 

2.2.1 Molecular Line Radiative Transfer 

We utilize the 3D non-LTE molecular line radiative transfer code, 
TURTLEBEAC H, to calculate the CO em ission properties from our 
model SMGs iNaravanan et alj|2008ah . TURTLEBEACH is an ex- 
act line transfer code that assumes f ull statistical equi l ibrium . The 
details of the code can be found in INaravanan et alj d2006bl) and 
INaravanan et alj d2008al) . and we refer the reader to these papers 
for details on the specific algorithms used. Here, we summarize the 
details relevant to this work. 

We assume that all cold, star-forming gas in the hydrodynamic 
simulations is neutral. Fol lowing the observational constraints of 
iBlitz & Rosolowskvl d2006h . we model the molecular fraction based 
on the ambient hydrostatic pressure: 



-Rmoi = nm/nm = 



3.5 x 10 4 



(1) 



where k is Boltzmann's constant, and the extern al pressure P ex t 
is calc ulated via the fitting formula derived by iRobertson et al.1 
d2004h : 

logP = 0.05(log n H ) 3 -0.246(log n H ) 2 +1.749(log n H )-10.6(2) 



d When considering emission properties which are strongly dependent on 
merger mass, we will, at times, employ the usage of SMG1 as a "high-mass" 
fiducial SMG along with SMG 10. This allows us to bracket the range of 
galaxy masses which form SMGs in our simulations. SMG1 forms an SMG 
ranging from average (S , g50~5-7 mJy during inspiral) to extremely lumi- 
nous and rare (Ss5o~15 mJy during final coalescence), whereas SMG10 
forms at peak (final coa lescence) only an average Ssso~5-l mJy SMG. See 
INaravanan et alj 120091) for more details regarding the lightcurve of model 
SMG1. 




Radius (kpc) 



Figure 2. Radially averaged molecular gas fraction as a function of galac- 
tocentric radius for three model SMGs (varying only the orientation angle 
of the disks) during three different phases. The three mod el SMGs produce 
"average" SMGs with fluxes ranging from S"850~2-7 mJy INaravanan et alj 
2009). The black lines are the pre-merger phase for each model; the blue 
lines denote the peak of the starburst (roughly the peak of the SMG phase); 
and the red lines represent the post-quasar phase. The molecular gas frac- 
tion is calculated via a dependency on the ambien t pressure as motiv ated by 
observations of clouds in local galaxies telitz & Rosolowskv 200(J). Dur- 
ing the peak of the SMG phase (blue lines), the bulk of the neutral gas in 
the central regions is molecular, owing to rather high gas densities in the 
galactic nucleus. 



In Figure |2j we show the radially averaged molecular gas fraction 
as a function of galactocentric distance for three model SMGsQ- 

In order to investigate numerous snapshots at relatively high 
temporal resolution (5 Myr) for a number of merger models, the 
hydrodynamic simulations were smoothed to a spatial resolution 
of ~ 160 pc. Resolution tests presented in lNaravanan et al] (T2008a) 
which investigated the spectral and spatial invariance of these meth- 
ods confirm that the resultant CO morphologies, lines profiles, and 
excitation conditions are convergent at these spatial resolutions. In 
order to properly account for the molecular density gradients which 
exist within cells of this volume (i.e. dense cloud cores and diffuse 
cloud atmospheres), we furt her model the molecular gas following 
the sub-grid prescriptions of Nar ayanan et all d2008al) . Specifically, 
the molecular gas within the ~160 pc grid cells is assumed to re- 
side in a mass spectrum of GMCs which are modeled as singular 
isothermal spheres. The mass spectrum of clouds follows a power- 



4 We note that thi s feature of the code is a n improvement over the algo- 
rithms described in Narayanan et al. 1 2008a). Specifically, previous works 
using TURTLEBEACH have assumed half the neutral gas (by m ass) was 
molecular, consistent with average conditions in local galaxies jKeres et"al] 
2003). While these assumptions tying the molecular gas fraction to global 
observations may r eproduce average molecular emissi on patterns in galax- 
ies and AGN (e.g. Narayanan et al. 2006a, 2008b. a c), it does not provide 
spatially resolved information regarding the molecular content. By tying the 
molecular gas fraction to the ambient pressure as motivated by observations 
of GMCs iBlitz & Rosolowsky|2006l) . we are able to more accurately model 
the spatial distribution of the molecular gas. We note that even more sophis- 
ticated models for treating the neutral gas breakdown ex ist in the literature 
(e.g. Pelupessy et al. 2006; Robertson & Kravtsov 2007). However, utiliz- 
ing methodologies such as these becomes computationally infeasible when 
considering the numbers of simulated galaxies and snapshots modeled in 
this work. 
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law with index a =1.8 as motivated by observations of local clouds 
jBlitz et alj2007l) , though tests have shown that TURTLEBEACH re- 
sults do not vary so long a s the mass spectrum indi ces reside within 
observational constraints dNaravanan e t al. 2008b). More details of 
the impl ementation of the s e sub-g rid treatments of GMCs may be 
found in iNaravanan et al. (2008a). Within the molecular gas, we 
conservatively assume that Galactic abundance patterns hold, and 
model the CO fractional abundance as CO/H2= 1.5 x 10~ 4 . While 
molecular abundances in the ISM of high redshift galaxies are 
unconstrained, recent measurements have shown that star-forming 
UV-selected galaxies, likely pro genitors of massive mergers, have 
solar abundances in their ISM (Shaplev et al. 2 0041) . As such, an 
assumption of Galactic molecular abundances may be reasonable. 

We build the emergent spectrum by integrating the equation of 
radiative transfer through the H2 gas: 

r 

I* = t 1 - ^ T ' Ar) ] e- T " (tot) (3) 

where /„ is the frequency-dependent intensity, S v is the source 
function, r is the physical distance along the line of sight, and r 
is the optical depth. 

The source function is made up of a combination of the emis- 
sion from dust, as well as line emission. Formally, S v = j v /a v 
where j v = j„(dust) + j„(gas). Similarly, a v has components 
from both the gas and dust. The dust radiates as a blackbody, and 
is assumed to be at the kinetic temperature of the molecular gas. 
We note that this element of the calculation is inconsistent with the 
SUNRISE calculations, which formally derives the dust temperature 
assu ming the dust and radiation fi eld are in radiative equilibrium. 
Here. lWeingartner & Drainej J200 lh opacities are assumed, though 
this makes little impact on the CO line flux. 

The line source function is dependent on the CO level popula- 
tions. Therefore, in order to self-consistently calculate the line in- 
tensities of CO, the molecular excitation properties must be known. 
The relevant physical processes in determining the CO excita- 
tion are collisions and radiative (de)excitation (e.g. line trapping; 
INaravanan e7ai]|2008ij) . The molecules are assumed to be in sta- 
tistical equilibrium, and the population levels are calculated con- 
sidering both the radiation field and collisional processes. 

The methodology is an iterative one. To determine the solution 
to the molecular excitation, the level populations across the galaxy 
are first guessed at (in practice, we guess a solution near LTE). The 
molecular gas is then allowed to radiate model photons based on 
the assumed level populations, and, when a sufficient number of 
photons have bee n realized in each grid cell, the mean intensity 
field is calculated dBernes|[l979l) . Under the assumption of statisti- 
cal equilibrium, the radiative excitation rates in combination with 
the collisional excitation rates give updated level populations, and 
new model photons are then emitted. This process is repeated until 
the level populations are converged. 

The line transfer takes into account velocity fields. The three 
dimensional velocity field across the model galaxy is returned by 
the GADGET- 3 hydrodynamic simulations. The line of sight veloc- 
ity gradient from cell to cell is accounted for in the line transfer via 
emission and absorption line profiles, both of which are Gaussian 
in nature. The emission and absorption profiles have their widths 
determined by the thermal line width in the cell, as well as an as- 
sumed microturbulent velocity field (here, set at 0.8 km s _1 ). The 
difference in the frequency centers of the emission and absorption 
profiles is determined by the line of sight velocity difference be- 



tween the emitting and absorbing clump (e.g. Equations 5 and 6 of 
INaravanan et alj2006bh . 

For the models presented here, ~13 million model photons 
were emitted per iteration. The mass spectrum of GMCs are con- 
sidered with a lower cutoff of lxlO 4 Mq, and upper limit of 
1 x 1O 6 M0 (con sistent with constraints provided by local GMCs; 
iBlitz et al]2 007). The CO excitation was solved for across 10 levels 
at a time, and the collisional rate coeffici ents were taken fro m the 
Leiden Atomic and Molecular Database dSchoier et al.ll2005l) . The 
boundary conditions included the cosmic microwave background 
which was modeled at 2=2.2 to have a temperature of T=8.74 K. 

Finally, we comment that because of spatial resolution lim- 
itations in the molecular line radiative transfer, we are forced to 
consider 8 kpc boxes. Rather than following the center of mass of 
the merging galaxy system (which, at times, may have scant little 
gas), we choose to follow a single progenitor galaxy through its 
evolution. Of course, as the nuclear disks of the progenitor galax- 
ies overlap (both during first passage, as well as coalescence), our 
models will include the emission from both galaxies in the pair. 

2.2.2 Polychromatic SED Radiative Transfer 

In order to identify when our simulated galaxies would be se- 
lected as submillimeter luminous sources, we simulate the ul- 
traviolet (UV) through submillimeter continuum photometry us- 
ing the 3D adaptive g rid polychromatic Monte Carlo radiativ e 
transfer code, SUNRISE d T onsson| [2006: Jon sson et ail 200(1 [2009). 
SUNRISE calculates the transfer of UV through millimeter wave ra- 
diation through the dusty interstellar medium. We refer the reader 
to these work for details on the underlying algorithms as well as 
numerical tests. Here, we summarize, and explain the physical pa- 
rameters employed in this study. 

The radiative transfer is implemented via a Monte Carlo al- 
gorithm in which photon packets representing many real polychro- 
matic photons propagate through the dusty interstellar medium, and 
undergo scattering, absorption, and reemission. Model cameras are 
placed around the simulated galaxy to sample a range of viewing 
angles. The emergent flux is determined by the number of photons 
that escape the galaxy unhindered in a given camera's direction, 
as well as those scattered into or reemitted by dust into the cam- 
era. For the purposes of the calculations presented here, we used 8 
cameras placed isotropically around the model galaxy. 

SUNRISE is able to handle arbitrary geometries for the sources 
and dust. The input spectrum includes contributions from stars 
and black holes (AGN). The AGN input spectrum util izes the 
luminosity-dependent templates of Hopk ins et al.l d2007l) of un- 
obscured quasars. The normalization of the input spectrum is 
set by the total bolometric luminosity of the central black hole, 
Lagn=??Mbhc 2 . Again, r/ is assumed to be 10% (c.f. § 12. It . 

The stellar input spectrum is calcul ated utilizing the stel- 
lar populations code, ST ARBURST99 dLeitherer et al] 1 19991 : 
Vazq uez" & Leitherej2005h . where the ages and metallicities of the 
stars are taken from the hydrodyna mic simulation s. We assume a 
Kroupa initial mass function (IMF) dKroupj|2002l) . consistent with 
recent resul t s from observations of z ~ 2 star form ing galaxies 
dDavell2008l : iTacconi et ai]|2008l : lvan"Dokkumll2008l) . The stellar 
particles initialized with the simulation are assumed to have formed 
over a constant star formation history. In order to match the star for- 
mation rate and stellar mass of the first snapshot of the simulation, 
this corresponds to a SFH of ~250 Myr. 

The stellar clusters with ages less than 10 Myr are assumed 
to reside in their nascent birthclouds. SUNRISE models the ef- 
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fects of reddening of the stellar spectrum through these birfh- 
clouds utilizing resu l ts from the photoionization code MAPPINGSIII 
dGroves et al.l |2004 l2008h . MAPPINGSIII calculates the transfer 
of continuum radiation and lines through the HII reg ions (which 
evolve as a one-dimensional mass-loss bubbles; ICastor et al.l 
and photodissociation regions (PDRs) surrounding stel- 
lar clusters dGroves et alj|2008h . The effect of modeling the ob- 
scuration of stellar clusters by HII regions and PDRs is a redis- 
tribution of emergent UV light into the far infrared/submillimeter 
bands. The HII regions absorb much of the ionizing UV flux and 
contribute heavily to the hydrogen line emission from the ISM as 
well as hot-dust emission. The time-averaged areal covering frac- 
tion by PDRs (/pdr) is related to the PDR clearing timescale as 
f pdr = exp(— t/tc iear), and is taken to be a free parameter (see 
iGroves et alj|2008l) . These PDRs absorb much of the non-ionizing 
UV radiation field, and contribute to the emergent PAH and FIR 
emission. Here, we assume an areal covering fraction of unity. A 
covering fraction ,f p dr=l translates to a clou d-clearing timescal e 
longer than the lifetimes of O and B stars dGroves et ai] |2008h . 
and the main consequence of reducing this clearing time scale is 
to red uce the emergent subm illimeter flux at the time of the star- 
burst dNaravanan et al.l2009T) . The cloud clearing time scales in gas 
rich galaxy mergers are unconstrained. That said, some constraints 
may be placed on this value either by educated guesses, as well 
as ansatzes that are then verified by a comparison of the simulated 
SEDs to those observed. 

The assumption of a clearing time scale longer than the life- 
times of O and B stars may be a reasonable guess. The cen- 
ters of local gas rich major mergers are known to have large 
molecular volume filling fractions, an d are well characterized by 
a uniform, smooth mo lecular medium dDownes & Solomonll 19981 ; 
ISakamoto etai1ll999h which may blanket nuclear O and B stars 
their entire lives. While an assumption of / p dr=l is not the same 
as a uniform molecular medium, tests have shown that in this lim- 
iting case, the submill imeter SED is the quite similar to / p dr=l 
dNaravanan et al.ll2009h . 

Similarly, we can take the reverse approach, and assume the 
ansatz of f pdr =l, and compare t he simulated SEDs to those ob- 
served. In lNaravanan et all d2009h . we showed that the mean SED 
in modeled SMGs with a range of masses compares quite well 
with observed SEDs of SMGs with spectroscopic redshifts. More- 
over, those authors found that the peak submillimeter flux was re- 
lated to the total mass of the galaxy. Utilizing PDR clearing time 
scales longer than the lifetimes of O and B stars resulted in average 
(S 8 50~ 5 mjy) SMGs with halo masses of order ~5x 10 12 - 10 13 
Mg, consistent with t he inferred halo masses of observed SMGs 
( Swin bank et ail l2008h . Similarly, the simulated stellar masses, 
black hole masses, and H2 masses both compared well with obser- 
vations, and scaled with peak submillimeter flux. Reducing t c iear 
(or, equivalently, / p dr) would require more massive galaxies to pro- 
duce relatively average Ssso ~ 5 mJy galaxies, and quickly violate 
the inferred halo masses from clustering measurements. As such, 
the assumption of / p d r = 1 may be reasonable. That said, it is im- 
portant to interpret the results presented in this paper in the context 
of this assumption for the birthclouds surrounding stellar clusters, 
especially in the context of mass-dependent results (such as the CO 
line widths; § !4.2t . 

The dust and radiation field are assumed to be in radiative 
equilib r ium , and utilize a method similar to that developed by 
lJuvelal d2005h . Here, when photons are absorbed in a grid cell, 
and the dust temperature updated, a new photon with SED equal 
to the difference between the SED emitted for the new dust tem- 



perature and that from the old one is emitted. This procedure is 
iterat ed upon until the radiation field has converged djonsson et al.l 
2009). The SUN RISE calculations em ploy 10 million photon pack- 
ets per iteration. I Jonsson et alj d2009f) find that the implementation 
of dust-temperature itera tion in SUNRISE recovers the solution to 
the lPascucci et alj d2004l) radiative transfer benchm arks to within a 
few pe rcent for UV through mm wavelengths. The lPascucci et all 
d2004h benchmarks were recovered through the most stringent test 
cases of T = 100. We note that the maximum optical depth we see 
during the SMG phase of these simulations is r ~ 75. 

The gas and stars initialized with the simulation are assigned 
a metallicity according to a closed-box model such that Z=(-y 
ln[/ gas ]) where Z is the metallicity, y the yield=0.02, and / gas 
is the initial gas fraction (though note the fluxes during the SMG 
phase of the model galaxies are not very sensitive to these assump- 
tions). Because the dust properties of z ~ 2 SMGs are relatively 
unconstrained, we utilize Galactic observations as input parame- 
ters. We assume a constant dust to gas ratio comparable to observa- 
tions of local ULIRGs of 1/50. Te ntative evidence suggests compa- 
rable dust to gas ratios in SMGs |Solomon & Vanden Boutl 1 20051 ; 
iGreve et alj[2005t iKovacs et al]|2006l ; iTacconi et aLlbOoiT simu- 
lations para meterizing the dust in terms of a Galactic dust to metals 
ratio of 0.4 dDwekll 19981) give similar results to within 10%. The 
dust grain model used is the R =3.1 Weingartner & Drainel d200lh 
dust model including updates bv lDraine&Lild2007t) . 

We are exploring the dust and photometric properties of z ~ 2 
galaxies in a number of companion papers at various stages of 
preparation. In this work, we focus on the CO properties of SMGs; 
as such, we primarily utilize SUNRISE as an informant as to when 
the model galaxies may have sub-mm luminosities comparable 
to observed SMGs. We classify our model galaxies as SMGs 
when they have an observer-frame 850/im flux ^ 5 mJy which is 
comparable to a 3- a detection in current wide-field surveys (e.g. 
Coppi n~e7al]|2006l) . We model our sources at z= 2.5, and there- 
fore this fiducial criterion corresponds to a flux limit at a rest-frame 
wavelength of A=243 /im. 



3 EVOLUTION OF SUB-MM, B-BAND AND H 2 
PROPERTIES OF SMGS 

We begin with a general description of the evolution of the 850 \im, 
B-band, and H2 properties of our model SMGs, and relate these to 
the evolutionary status of the galaxy merger. In Figure[3] we present 
the 850 p_m flux, CO line width (which will be discussed more in 
§ 14. 2\ , and observed B-band apparent mag nitude (AB magnitude 
system , modeled at £=2.5, typical of SMGs; [Chapman et al.l2003al . 
l2005h of our fiducial merger simulation, SMG10 (Table[2j}- The 
blue shaded region represents the sightline dependent range of po- 
tential B-band magnitudes. The CO transition plotted is J=3-2 in 
the rest-frame, corresponding to mm-wave observations at z ~ 2. 
As a reference for the global morphology, in Figure [4] we plot 
the projected gas density at various snapshots for fiducial model 
SMG 10, with the location of the black holes overlaid. 

The initial passage of the galaxies induces a starburst of order 
~200 MQyr -1 . The galaxies form stars at this rate for a few x 10 8 
yr as they inspiral toward final coalescence. Duri ng this time, the 
galax y builds a stellar mass of order ~10 11 Mq dNaravanan et al.1 
2009). This elevated SFR allows the galaxy to produce ~l-2 mJy 
of flux at 850 p.m rendering this galaxy below the nominal 5 mJy 
detection threshold for SMGs, though detectable with future sensi- 
tive instruments. 
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Figure 3. Evolution of 850 ytm flux (top row), CO (J=3-2) line width (second row) and observed-frame B-band flux (AB magnitude; third row) for model 
galaxy SMG10 (Table[2j] designed to be an average SMG). All quantities are plotted at redshift z=2.5. The grey hatched region in the top plot shows when 
the galaxy is visible as an SMG above Ss50> 5 mJy. The red shaded region in the CO line width plot shows the dispersion over 100 random sightlines, 
and the blue shaded region in the B-band plot denotes the dispersion over 8 cameras placed isotropically around the galaxy. As the galaxies spiral in toward 
coalescence, the observed submillimeter flux is relatively low (S85o~l mJy), and the progenitor galaxies are disk-like. Consequently, the line widths are 
representative of the virial velocity of the galaxies (here, set at V^=320 km s _1 ). As the galaxies coalesce (T 0.6-0.65 Gyr), the ~1000 MQyr -1 starburst 
drives the 850 (am flux to detectable levels (/Sg5o> 5 mJy). Concomitantly, the B-band flux rises sharply, both owing to the intense starburst and a rapidly 
growing AGN. The CO line FWHM doubles as two disks enter the simulation box/observational beam (here, set at 8 kpc). In the post-starburst phase, the 
submillimeter and B-band flux drops, and the line widths settle toward the rotational velocity of the combined (two-galaxy) system. See main text for more 
details. 
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Figure 4. Projected gas density for a fiducial merger model SMG10 (designed to be an average SMG). Box sizes are 100 kpc on a side, and the box centers 
on the same galaxy the CO calculations center on. The time stamps are in the upper left of each panel, and in units of Gyr. For reference, the SMG phase is at 
time T 0.65 Gyr. Time stamps near the SMG phase have panels with equivalent time stamps to Figures|5]and [7] For these panels, the yellow box highlights 
the 8 kpc simulation box employed for the CO radiative transfer calculations. The location of black holes is noted by the red crosses. 
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When the galaxies approach toward final coalescence (T~0.6 
Gyr), tidal torques from the merger drive large-scale inflow s 
dBarnes & HemquisJl99llll996l : lMihos & Hernquistlll994l . ll996h . 
Physically, the interaction of the galaxies spin up the disks as 
they transfer angular momentum from the orbit of the galaxies 
to the disks themselves, and triggers the growth of bars. Gas that 
shocks on this bar dissipates energy and loses angul ar momentum, 
causing an inflow to the central regions (see, e.g. Ho pkins et al.l 
l2008t 2009). The high densities in the nuclear region of the co- 
alesced system give rise to a massive starburst of order ~1000- 
1300 MQyr -1 . The nascent PDRs surrounding young stellar clus- 
ters convert the UV flux intercepted from O and B stars into longer 
wavelength submillimeter radiation. During this starburst event, the 
galaxy may be selected as a luminou s submillimeter source with 
850 |i.m fluxes approaching ~5-7 mJv ( Nar ayanan et al.l2009t> . The 
peak submillimeter flux observed is directly related to the mass 
of the merging galaxies - galaxy mergers above total (halo) mass 
~5x 10 12 Mq will produce the nominal ~5 mjy at 850 \im to be 
detectable as an SMG. Mergers of a significantly lower mass will 
have difficulty producing a strong enough starbur st to drive the ob- 
served submillimeter flux (Naray anan et alj20 09). This owes to the 
fact that the submillimeter flux in our model derives largely from 
reprocessing of UV flux from the starburst i n the birth clouds sur- 
rounding stellar clusters ^Groves et al . 2008). 

Concomitant to the final coalescence starburst, inflows fuel 
central black hole accretion. A fraction of the accreted mass energy 
is deposited into the ISM surrounding the central black hole(s), 
driving a pressure-driven wind. These AGN winds expel much of 
the obscuring gas and dust in the central regions, allowing the sys- 
tem to be viewed as an optically selected quasar ( T~0, 65-0.7 Gyr; 
e.g. lHopkins et al]2005l . [2006l ; ISpringel et al.l2005l , and references 
therein). The quasar phase and SMG phase are roughly coincident 
(thoug h the quasar phase may lag the SMG phase by up to ~20 
Myr; Springel et al. 2005). This is in good agreement with obser- 
vations bv lCoppin et al.l d2008bl) , who find an overlap in a subset of 
their z ~ 2 observed SMGs and quasars. As demonstrated by Fig- 
ure [3] however, there is a large sightline-dependent dispersion in 
potential B-band fluxes during the final merger. Consequently, the 
same galaxy can show almost 2 magnitudes dispersion based on 
observed viewing angle, and not all SMGs will appear as quasars. 
However, the galaxy is visible as an SMG at all modeled sightlines. 

It is important to note that the selectability of our simu- 
lated galaxies as quas ars is mass-dependent. As discussed by 
iNaravanan et alj ^20091) , the final black hole mass of the merged 
system is dependent on the m a ss of the g alaxy. Th i s is s imilar 
to results found by llidz et al.l fcOOd) and lli et al.l fe007l) who 
found that quasar lumin osity is tied to progenitor galaxy halo mass. 
INaravanan etall J2009h found that the final black hole masses 
of average SMGs (e.g. S850«5-7 mjy) will be a fewxlO 8 M Q , 
whereas the mergers which produce the most luminous SMGs 
(<S l 850~15-20 mjy) may make black holes comparable to those seen 
in quasars (A/bh ~10 9 Mq). We therefore continue the discussion 
of the CO properties of our SMGs during the "quasar phase" as 
referring to the time period when the simulated dust-attenuated B- 
band flux peaks (e.g. third panel, Figure(3](, and note that this may 
not necessarily correspond to a galaxy selectable as a quasar in cur- 
rent surveys. The exact relationship between SMGs and quasars 
will be discussed in due course (J. Younger et al. in prep; D. 
Narayanan et al. in prep). 

The starbursts induced by the merger consume significant 
amounts of gas. However, while gas consumption is high, su- 
pernova pressurization of the ISM sustains large molecular gas 
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Figure 6. Relative fraction of gas in rotational motion about the galactic nu- 
cleus (solid line, units on left axis), and 850 |jm flux (dotted line, units on 
right axis) for fiducial merger model SMG10. Following Cox et al. (2009, 
submitted), gas is considered to be in in rotational motion when it has cir- 
cular velocity at least 75% of the expected Keplerian velocity at its radius. 
During the peak of the SMG phase, most of the gas is disturbed from the 
central rotationally supported disk. 



fractions. Thermal energy input into the hot phase ISM, as well 
as mass increases owing to the stellar mass returned and evap- 
oration of cold clouds increases t he ambient pressure on cold 
clouds JSpringel & Hernquistll2003l) . This increase in pressure in- 
creases molecular fractions in the neutral ISM, in accordance with 
our pressure-based H2 formation/destruction algorithm (c.f. Equa- 
tion 0. Hence, while star formation of course consumes H2 gas, 
during starbursts, this effect may be mitigated owing to conversion 
of HI to H2. The H2 masses during the final coalescence burst are 
rather high with typic al masses ^5x 1 10 M n , though can be as 
high as ~3x lO n M ( INaravanan e"tai]|2009l) . 



4 THE EVOLUTION OF CO PROPERTIES OF SMGS 
4.1 CO Morphologies 

4.1.1 Molecular Disk Formation and Disruption 

The structure of the molecular gas in SMGs, and in particular the 
(potential) existence of molecular disks is essential to a thorough 
understanding of the evolution of their CO line properties. We 
briefly outline the key points related to this topic here, though defer 
a detailed investigation into the survivability of molecular disks in 
mergers to a future work. 

In Figure [5] we show the CO (J=3-2) morphology of galaxy 
SMG 10 as a function of time. The temporal evolution depicts the 
CO morphology as it evolves through inspiral (first row) and final 
coalescence (second and third row). The galaxy is most likely to 
be viewed as an SMG during fi nal coalescence when t he SFRs are 
most elevated (e.g. Figure 1, INaravanan et al.l [2009T) . Here, this 
corresponds to T~0.65 Gyr. Because our CO calculations center 
around a single galaxy at all points in its evolution to maximize 
spatial resolution, both galaxies only appear in the images when 
the nuclei are both within the 8 kpc model box (second row and 
beyond). 

In Figure [6] we show the relative fraction of gas in rotational 
motion as a function of time, with the sub-mm flux overlaid as a ref- 
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Figure 5. Simulated CO (J=3-2) morphology for fiducial model SMG10. The boxes are 8 kpc on a side, and the intensity is in velocity-integrated brightness 
temperature with scale on bottom. The simulations focus on a single galaxy through its evolution until both galaxies are within a single 8 kpc box (second row 
and beyond). For reference, the 8 kpc box employed for the CO radiative transfer simulations is shown explicitly with respect to the global morphology in 
Figure|4] During inspiral (first row of this figure) the molecular gas is in a disk-like configuration. As the second galaxy enters the box, and they merge (second 
row), the disk-like morphology is disturbed and a large fraction of the gas is pulled from disk-like motion into relatively radial orbits. This corresponds to the 
peak of the SMG phase. During this time, extended features and tidal tails may become apparent (middle row). Tidal torquing drives much of the gas toward 
the central regions, resulting in a relatively concentrated molecular gas spatial extent (third row). 
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erence. Gas is considered to be in rotational motion when its circu- 
lar velocity exceeds 75% of that expected at its radius for Keplerian 
motiorQ. As such, the relative fraction of gas in rotational motion 
may be viewed as a measure of the "diskiness" of the molecular 
gas. Comparisons between Figures [5] and [6] may be made via the 
time stamps displayed in the panels of Figure [5] 

The galaxies remain relatively disk-like as they inspiral toward 
final coalescence. Upon final merging, when the galaxy undergoes 
its luminous SMG phase, the disk is tidally disturbed. Tidal features 
(and, on occasion, double-nuclei) become apparent in the molecu- 
lar gas morphology during this final-coalescence SMG phase (e.g. 
T~ 0.65 Gyr, Figure [5}. This can be seen more explicitly in Fig- 
ure[7] where we plot the CO (J=3-2) centroid velocity maps of the 
same snapshots shown in Figure [5] Soon after the final interac- 
tion/SMG phase, the gas yet again re-virializes, and a strong (com- 
pact) molecular disk re-forms (this is seen in Figure [6] though the 
snapshots in Figure [7] do not extend far enough in time to show 
this phase). This history of molecular disk formation/disruption 
throughout the galaxy merger's history is reminis cent of that seen 
in mod els of the molecular ISM in z ~ 6 quasars dNaravanan et al.1 
2008c), and will play an important role in our understanding of the 
CO line widths, and usage of CO as a dynamical mass indicator in 
the sections to come. 



4.1.2 Spatial Extent of CO Emission 

With recent advances in (sub)mm instrumentation, high spatial res- 
olution images of t he molecular gas in SMGs is beginning to be- 
come available (e.g. lTacconi et alj2008l) . Constraints on the spatial 
extent of the CO emission from SMGs can be important for e.g. 
determining the dynamical mass from CO line widths. 

During the merger-induced starburst, tidal torquing drives cold 
gas from the outer disk into the nuclear regions of the galaxy. 
Consequently, the CO emission becomes rather compact. In Fig- 
ure [8] we plot the distribution of characteristic CO (J=3-2) radii 
from model galaxies SMG1 and SMG10 over 100 random sight- 
lines during snapshots where the system may be viewed as an SMG 
(£850 > 5 mjy). The characteristic radius for CO emission is the 
flux- weighted standard deviation in the radial distribution of fluxes. 
Practically, the map is treated as a histogram of fluxes at varying 
radii. The standard deviation in this histogram is the characteristic 
radius. 

During the SMG phase, the emission is relatively compact ow- 
ing to the gas funneled into the central regions from the merger. 
The characteristic radius averages at ~ 1.5 kpc for most sightlines, 
with some dispersion owing to evolutionary status and sightline. 
Some large CO radii are seen from inspiralling disks in massive 
(Mum ~ 10 13 M,-^) mergers which a re already SMGs during the 
inspiral phase dNaravan an et al. 2 0091) . 

The CO spa tial extents modeled here are comparable to the 
measurements bv lTacconi et al .1 d200q) of a ~4 kpc FWHM diame- 
ter in SMGs (which corresponds to a ~0.85 kpc standard deviation 
in radius if the emission is Gaussian in nature). The distribution 
of radii additionally is consistent with the a ~ 0.6-3.3 kpc range 
of spatial extents observed from SMGs dDownes"& Solomo nl2003l ; 



5 Varying this fiducial fraction of 75% makes no difference on the relative 
temporal evolution of the fraction of disk-like gas; lowering or increasing 
this value simply increases or lowers the normalization of the curve. As 
such, throughout this paper, this "fraction" should be taken as relative, and 
not absolute. 
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Figure 8. Characteristic radius of CO (J=3-2) emission during the SMG 
phase of fiducial models SMG1 and SMG 10. These models were cho- 
sen to bracket the range of masses which p roduce SMGs ( ranging from 
average, (5-7 mjy), to luminous (~15 mjy); [Narayanan et al. 200^). The 
characteristic radius is calculated using the standard deviation in the flux- 
weighted radial distribution of fluxes. The distribution is plotted over 100 
random sightlines. The charact eristic radii match up well with observed 
size scale FWHM measured by Downes & Solomon! i2003); Genz el et alj 
l2003l) : lTacconi et alj <2006ll2008h and lYounger et al.1 J2008I) . Please note, 
however, that because the simulations here are performed in isolation, rather 
than drawn from cosmological conditions, this distribution represents the 
range of expected CO (J=3-2) spatial extents from observed SMGs, rather 
than a true distribution. 



Genzel et alj|2003l : iTacconi et af1l2006l . 120081; lYounger et ai]|2008l ; 
Iono et al.ll2009l) 



4.2 CO Line Widths 

4.2.1 Model Results: Evolution of CO Line Widths 

The observed molecular line widths of SMGs are exceptionally 
broad, with a median FWH M of ^600-800 km s -1 , and line widths 
exceeding 1000 km s" 1 dGreve et alJl2005L ICarilli & Wandl200rj; 
lTacconietal]|2006l ; ICoppin et alj|2008bl : llono et alj|2009l) Here, 
we explore the evolution of CO line widths in our model SMGs; 
we show how our merger-driven formalism for SMG formation and 
evolution may self-consistently explain the observed broad lines 
from SMGs. Because the CO emission line is essentially a distribu- 
tion measuring the power at a range of molecular gas line of sight 
velocities, rather than employing any particular fitting methodol- 
ogy (and thus espousing the associated uncertainties), we treat the 
line as a distribution of fluxes, and utilize the standard deviation of 
the distribution (a) as a measure of the line width. For a perfectly 
Gaussian line, the FWHM of the line would simply be ~2.35 xer. 
While the lines are not perfectly Gaussian, to better compare with 
observations, we utilize this conversion between a and FWHM. 

In short, the CO line widths are repre sentative of the dynam- 
ics of the system dNaravanan et al.|[2008cl) . Prior to final coales- 
cence, typically only a single galaxy is within the beam, and the 
line widths are narrow, representing the rotational velocity of a sin- 
gle galaxy. During the SMG phase, when the galaxies merge, the 
line widths roughly double owing to the contribution of emission 
from both disks. In the post-SMG phase, as the molecular gas re- 
laxes into a new disk (Figure [5), the line widths drop by a factor 
V2. 
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Figure 7. Simulated CO (J=3-2) centroid velocity maps for model SMG10. The boxes are 8 kpc on a side, and the velocity are in units of km s with scale 
on bottom. The snapshots match those of Figure [5]The inspiral phase (top row) is characterized by ordered disk-like motion. The molecular disks are rapidly 
destroyed as the galaxies coalesce during the SMG phase (T = 0.63 — 0.65 Gyr). Note, the disk angle slowly changes throughout the early part of the galaxy's 
evolution, thus changing the magnitude of the line of sight velocities seen. 
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In more detail, the CO line widths reflect the dynamics of the 
molecular gas in the galaxy. This owes to the origin of the CO emis- 
sion lines from the model galaxy. While emission within a given 
cell (or neighboring cells) is typically optically thick, globally, the 
emission is optically thin. When the CO line emitted from a cell 
escapes the local region, it typically leaves the galaxy unhindered 
as steep velocity gradients in the throes of the merger shift the ab- 
sorption profile out of resonance with the emission profile. The CO 
emission line, then, is essentially the sum of the emission that es- 
capes individual cells containing GMCs at their individual veloci- 
ties. As such, the line widths are reflective of the dynamics of the 
system. The dispersion (<r) in the CO line width corresponds well 
with the velocity dispersion of the gas along the line of sight. 

As the galaxies inspiral toward final coalescence, the molecu- 
lar gas in the progenitor galaxies of the SMG is relatively virialized 
(Cox et al., 2009, submitted), and thus the line widths are reflec- 
tive of the circ ular velocity of the pro genitor disks. As discussed 
in Table |2. ll and lNaravanan et all ( 2009). our model mergers which 
produce average SMGs (S850~5-7 mjy) were typically initialized 
with disks with K=320 km s _1 , thus resulting in a CO linewidth 
of a w 160 km s _1 (which is equivalent to the 320 km s _1 circular 
velocity in the disk modulated by sin(30° ) to account for the aver- 
age inclination angle of the disk). This corresponds to a FWHM of 
~375 km s~ for a Gaussian velocity dispersion (Figure[3j. Recall 
that during the inspiral phase, typically only one galaxy is within 
the 8 kpc beam, which is comparable to mos t interferometric beam 
sizes jGreveetal.f2 005; Tac coni et al]|2006h . 

When the galaxies coalesce, both galaxies contribute to the 
detected line. The velocity dispersion of gas detected doubles, and 
consequently, the line width roughly doubles. This is coincident 
with the SMG phase. Here, this causes the line widths to increase 
to 2 x V c . Because the FWHM « 2.35 x V c , and accounting for the 
average inclination angle of the disks, we then arrive at a general- 
ized expression for the CO line FWHM during the final coalescence 
SMG phase: 



FWHMsmg ~ 2 x sin(i) x 2.35 x V c 



(4) 



which is, of course, simply FWHM w 2.35 x V c when an average 
inclination angle of i=30° is assumed. Here, where circular veloc- 
ities of V c =320 km s~ are employed for the initialization of the 
progenitor galaxies, this results in modeled linewidths of ~600- 
800 km s~ . The dispersion in the modeled CO line widths owes 
to both viewing angle and evolutionary effects. 

Generically, the SMG phase and quasar phase occur at or 
around the nuclear coalescence of the galaxies. However, the exact 
timing of this event is dependent on a number of things, includ- 
ing galaxy orbit, gas content, and mass, amongst other factors. As 
such, the CO line widths from the SMGs or quasar host galaxies 
with a merger-origin show a large range of line widths. In Figure[9] 
we plot the histogram of CO linewidths from our fiducial galaxy 
(SMG10) during both its SMG phase as well as its quasar phase. 
As before, we nominally define the SMG phase as when the 850 
p.m flux is > 5 mjy, and arbitrarily define the quasar phase as when 
the galaxy is brighter than 27th magnitude (AB Magnitude; appar- 
ent magnitude at z=2.5). Because the SMG phase and quasar phase 
are nearly coincident, the linewidths are generally broad throughout 
both phases, with nearly indistinguishable distributions. 

4.2.2 Observational Comparisons 

The distribution of modeled line widths for SMGs shown in 
Figure [9] is in excellent agreement with those presented by 



ICoppin et al] J2008bl) and ICarilli & Wand <2006h . The linewidths 
are broadly reflective of the mass of the simulated SMGs, and likely 
signifies a correspondence between the final masses and evolution- 
ary status of our modeled SMGs and those in nature. In this sense, 
our model faithfully provides a natural explanation for the broad 
line widths observed in SMGs. 

A comparison between the line w idths of our model q uasars 
and those in nature is more difficult. ICoppin et al.l d2008bl) pre- 
sented new CO detections of submillimeter-luminous quasars, and 
found a nearly indistinguishable distribution of CO line widths 
from quasars and SMGs, in excellent agreement with the simu- 
lat ions presente d here. In both our models, and the observations 
of ICoppin etal] d2008bh , the median CO line width is ~600-800 
km s _1 , which we view as a general success of our model. How- 
ever, a lit erature compilat i on of line widths from CO-detected 
quasars bv lCarilli & Wand J2006h found a much narrower median 
in the distribution, with the line widths showing a median value of 
F WHM~300 km s -1 , in contrast to both the observational results 
o flCoppin et all (2008b ) as well as the model results presented here. 
The source of this is discrepancy not entire ly clear. Because the 
analysis performed bv lCarill i & Wans ( 2006) was from a literature 
compilation, the sample utilized a variety of CO transitions from 
a non-uniform sample with objects spanning a large range of red- 
shifts from 2 ~ 1-6. Different CO lines trace different spatial ex- 
tents, and thus may show a range of line wid ths. Similarly, quasars 
even of a similar mass as those presented in Coppin et al. (2008b) 
but at a lower redshift would naturally show evolution in their line 
widths owing to shallower gravitational potentials. In contrast, the 
observations by Coppin et al. utilized only CO (J=2-l) and (J=3- 
2) emission from a sample of quasars in a smaller redshift range 
(«=1.7-2.6). In this sense, the comparison between our models and 
the relatively uniform observations of z ~ 2 quasars by Coppin et 
al. seems most appropriate as it best compares to the redshifts and 
emission lines investigated in this work. 



THE USAGE OF CO AS A DYNAMICAL MASS 
TRACER IN SMGS 



In § 14.21 we saw that during the SMG phase of the model galaxy's 
evolution, the CO line widths were larger than expected from or- 
dered disk-like motion gas. In this context, it is interesting to quan- 
tify the usage of CO as a dynamical mass tracer in SMGs. 

The exact value of CO line widths in high redshift mergers as 
a dynamical mass indicator depends, of course, on the relationship 
between linewidth and dynamical mass assumed. Broadly, for disk- 
like motion, this relationship can be expressed in the form: 



-^dyn — h 



^^fwhm^liwhm 

Gsin 2 (i) 



(5) 



where a is the ID velocity dispersion in the line, i the disk inclina- 
tion angle, R the spatial extent of the CO emission (here, taken to 
be the half-width at half-maximum), and k a constant encompass- 
ing the relationship between V c and Vfwtm (for a given distribution 
of mass) and R g and i?hwhm- Equation [6] holds due to the global 
optical thinness of molecular line emission across a galaxy. 

However, uncertainties exist in all of these conversion factors. 
Furthermore, it is not clear what fraction of the molecular disk sur- 
vives during the merger, and whether the inclusion of a sin 2 (i) term 
is appropriate. Even if it were, prior to the high resolution imag- 
ing capable only by ALMA, we can at best assume an average disk 



14 Narayanan, Cox, Hayward, Younger & Hernquist 



o 
z 



E 




angle of i=30°. As such, a constructive method for quantifying the 
usage of CO as a dynamical mass tracer in SMGs is to say: 
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FWHM (km s~') 

Figure 9. Distribution of modeled CO (J=3-2) linewidths during the SMG 
and "quasar" phases of model SMG 10. The SMG phase is considered when 
>S I 850>5 mJy, and the quasar phase is arbitrarily assigned with an appar- 
ent 73-band magnitude (AB system) cut of 27th magnitude. See Figure [5] 
for more details. The phases occur at similar time periods (typically sepa- 
rated by at most ~20 Myr), and consequently have similar line width dis- 
tributions. The FWHM distribution for these two sources appear to cor - 
respond well with the distributions published by Coppin et al. 1 2008b), 
th ough the quasar d istribution appears to be discrepant with those published 
bv lCarilli&Wand feOOd) . See text for details regarding the origin of the 
broad CO line widths in SMGs and quasars. Again, because the simula- 
tions here are not cosmological in nature, the distribution should be viewed 
as representative of the range of CO line widths from average SMGs, not 
the true distribution from a blind survey. While we normally include higher 
mass models (e.g. SMG1) in our analysis, we refrain here as the extremely 
broad (FWHMpsIOOO-1400 km s" 1 ) lines from SMG1 take away from the 
main point of modeled SMGs with average Ssso fluxes producing average 
CO line widths. SMGs as massive as SMG1 will be rare (indeed, this is ap- 
parent by the relative lack of observed Ssso > 15 mJy SMGs), though when 
observed, will have line widths FWHM > 1000 km s — 1 on average. 




Figure 10. Sightline averaged distribution of M dyn /M actU ai for models 
SMG1 and SMG10. These models were chosen to bracket the range of 
masses of galaxies which appear to produce SMGs in our simulations. The 
dynamical masses are calculated as in Equation [5] for the simulated CO 
(J=3-2) emission. Generally, using CO line widths from SMGs for a dy- 
namical mass calculation requires a modulation of a factor of ~ 1.5-2 for 
translating to a true enclosed mass. See text for details. 



M, 



dyn — 



^^fwhm^hwhm 



G 



= k' x M a 



(6) 



and characterize the relationship between AV f ^, hul Rh w hm/G and 
M ac tuai. Here, note that M ac tuai is the total (baryonic and dark 
matter) mass enclosed in 7?hwhm, and AVf w h m is calculated as 
2.35 x the standard deviation in the line width (c.f. § |4.2j . In Fig- 
ure [10] we plot the ratio of the dynamical mass as calculated in 
Equation [6] to the true enclosed mass (which includes dark mat- 
ter and baryonic mass). The radius employed is the half-width at 
half maximum from the simulated CO emission maps. The ratio 
of Mdyn/Mactuai is calculated for SMGs which span the range of 
masses explored here, for all snapshots which satisfy the fiducial 
criteria Sg5o>5 mJy, and for 100 random sightlines. On average, 
the ratio between the dynamical mass inferred from Equation 6]and 
the true enclosed mass in the half-mass radius, k', range^j from 
1.5-1.9. With this, it is then feasible to modulate Equation [6] with 
this given value of k' to construct the appropriate relationship be- 
tween CO-derived dynamical mass and true mass enclosed in the 
CO emitting region. 



6 CO EXCITATION AND LINE SPECTRAL ENERGY 
DISTRIBUTION 

The excitation of CO is an important diagnostic of high redshift 
galaxies. First, molecular line measurements from SMGs are typ- 
ically made in millimeter bands, which corresponds to high exci- 
tation CO lines in the rest frame at z ~ 2. Because the inferred 
molecular gas mass is generally derived via converting CO (J=l- 
0) velocity integrated intensity, understanding the typical CO line 
ratios (and average level of thermalization of the gas mass within 
the telescope beam) is crucial. Assumptions of thermalized line ra- 
tios (e.g. LTE, when brightness temperature ratios between levels 
are unity) between higher lying lines and CO (J=l-0) may under- 
estimate the molecular gas mass in the event of substantial quan- 
tities of subthermal gas. Second, the CO excitation patterns re- 
veal the line(s) of dominant CO power output. As more broadband 
molecular line spectrometers beco me available (e.g. the ZEUS, Z - 
spec, Zpectrometer spectrometers; [Hailev-Dunsheath et alj d2008l) : 
iBradford et"ai] ( 120041) ; lHarris et al] d2007l) ). CO detections of high 
redshift galaxies will be critically dependent on observations of the 
brightest lines. In an effort to aid interpretation of existing data sets, 
and help guide future observations of SMGs, we investigate the 
molecular excitation properties and CO line ratios from our model 
SMGs. We do not attempt to address the long standing problem of 
converting CO (J=l-0) flux to H2 gas mass in galaxies, but rather 
simply relate the flux density from higher excitation lines to that 
from the ground rotational transition. 



6.1 Model Results: Highly Excited CO in SMGs 

In Figure QT| we plot the sightline-averaged model CO SED for 
models SMG1 and SMG10 (Table Hj} which bracket the mass 
range of galaxies that satisfy the nominal selection criteria 5*850 > 



' Because these simulations are not cosmological, it is not possible to quote 
an exact mean in the distribution of AfdynAMactual . This range signifies the 
range seen when varying the exact model SMG or number of SMGs used 
to create a distribution similar to that plotted in Figure ITol 
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Figure 11. Model CO SED from fiducial models SMG1 and SMG10. The SED is presented for all snapshots where the galaxy would be detected as an 
SMG (5850 > 5 mJy), and the solid line shows the mean SED while the grey shaded region denotes the 1-cr dispersion amongst snapshots. The flux densities 
from J>1 levels are compared against the CO (J=l-0) level. The dashed line represen ts the pre dicted CO SED for fhermalized level populations (LTE). 
The g reen diamond show s observational data from SMG SMM 13120+42 42 Mainline et a?]|2006l) , and the red triangles from SMG ERO J 164502 (aka HR 
10) jAndreani et al.ll2000l;lGreve et alJl2003l;|Papadopoulos & Ivisonll2002h . To our knowledge, these represent the only SMGs with tabulated multi-line data 
(including a CO (J=l-0) detection) which are unlensed. The highest excitation CO SEDs come from final coalescence mergers, whereas massive SMGs that 
are above 5g5o>5 mJy during inspiral may have lower CO excitation. The mean SED displays reasonable agreement with the observations, and suggest that 
CO levels J > 2 may not be fhermalized. This shows that caution must be exercised when assuming a brightness temperature ratio of unity to the CO (J=l-0) 
line when deriving molecular gas masses. 



5 mJy. The ordinate is the ratio of the velocity-integrated intensity 
from various transitions compared to CO (J=l-0), normalized to 
CO (J=l-0). The range of emergent CO SEDs from our models is 
denoted in the grey shaded region, and the mean by the solid line. 
The dashed line shows the expected CO SED for thermalized level 
populations. The line ratios are modeled for unresolved observa- 
tions, and include emission from the entire 8 kpc simulation box. 
To best compare with the few published constraints (next section), 
we plot the line ratios relative to the ground (J=l-0) transition. This 
has the added benefit of providing a direct measure of the ability of 
higher lying transitions to convert to H2 gas masses. 

As Figure QTJ shows, there is a broad dispersion in potential 
CO SEDs from our simulated SMG. However, two generic features 
are evident. First, on average, the CO is quite excited with the SED 
turnover only occurring at the J~5-6 level. This is in reasonable 
agreement with observations of z ~ 2 SMGs (see next section). 
The higher excitation CO SEDs come from SMGs which arise dur- 
ing final coalescence mergers, while the lower excitation SEDs rep- 
resent the inspiral phase of massive galaxies (e.g. SMG1) which 



m ay be seen as SMGs ev en prior to final coalescence (e.g. Figure 
1, iNaravanan et al.ll2009l) . 

Second, while the turnover is typically at a relatively high J 
level, most CO levels are subthermal over the 8 kpc simulated box. 
The rotational ladder in Figure QTJ show that the gas is nearly ther- 
malized for the lowest lying lines (CO J=2-l). Higher lying lines, 
however, are not thermalized over the 8 kpc model box. While the 
gas associated with the nuclear starburst in the central ~1 kpc is 
warm, dense, and thermalized, the outer regions (R~l-4 kpc) con- 
tain a significant quantity of lower density gas which contributes 
to the emergent spectrum via line trapping. This lowers the mean 
excitation conditions observed. As such, Figure QTJ demonstrates 
that most observations of SMGs which do not focus solely on the 
nucleus will reveal subthermal CO emission at higher-lying (J>3) 
levelfl This is effectively saying the filling factor of dense gas is 
relatively low, thus lowering the mean observed CO SED. Conse- 



Higher spatial resolution observations probing the nuclei of SMGs should 



16 Narayanan, Cox, Hayward, Younger & Hernquist 



quently, assumptions of thermalized CO line ratios from unresolved 
observations ofSMGs will typically underpredict the molecular gas 
mass. 

6.2 Observational Comparisons 

While few tabulated observational constraints exist for relative in- 
tensities from SMGs, the excitation seen in Figure QT] seems to 
display reasonable agreement with the existing published data. 
Line ratios for ERO J 164502 and ERO J 164502 a re given by 
lAndreani et alj 12 000); Papa dopoulos & Ivisonl d2002h : lGreve et all 
j2003l) and lHainline et all t2006). These data are represented by the 
points with error bars in Figure [TT] While the CO SED is normal- 
ized to the ground state (thus making this point an unconstraining 
comparison), the higher excitation detection of two lines (J=2-l and 
J=4-3) is entirely consistent with the predictions made here. The 
intensity from the J=5-4 line from ERO J164502 is lower than the 
models predi ct here, though we note that this is lowest exc itation 
SMG known dPapadopoulos & Ivisonll2002l ; IWeifi et all2007l) . 

A detai led survey and com pilation of literature data has been 
presented bv lWeiB et alj d2007h . While these data are not tabulated, 
we can make qualitative comparisons. The bulk of these galaxies 
show subthermal emission in the higher lying lines with a turnover 
at the J=5 or 6 level. Thus, at face value, these observed CO SEDs 
are comparable with the predictions made in Figure [TT1 

A detailed examination of the WeiB et al. observed CO SEDs, 
however, suggests a potential discrepancy between the observed 
values and those modeled here. The observed intensity at the J=5 
or 6 level (where the SED turns over) is a factor of ~15 above that 
from the ground (J=l-0) transition; this is in contrast to our models 
which suggest that the most excited lines will be only a factor of 
~5 greater than the ground state (Figure [TTt . This difference may 
be reconciled with better constraints on the true CO (J=l-0) emis- 
sion from SMGs. Recall that the CO SED is typically plotted (both 
here, as well as in observational literature) as relative to the CO 
(J=l-0) transition. For the bulk of the observed sources in the liter- 
ature, the J=l-0 line has not directly observed, but rather inferred 
from large velocity gradient (LVG) modeling which is simplified 
approximate method compared to the 3D non-LTE radiative trans- 
fer methods employed here. It may be that the apparent difference 
in excitation between the observations and theoretical predictions 
will be reconciled via future detections of CO (J=l-0) from SMGs. 



7 DISCUSSION: TESTABLE PREDICTIONS 
7.1 Testable Predictions 

In this paper, we have outlined a model for the CO emission from 
high redshift SMGs. Our models have provided a natural explana- 
tion for the observed H2 gas masses, CO spatial extents, line widths 
and excitation conditions. We have attempted to compare with lit- 
erature data when available, and generally found a reasonable cor- 
respondence between our model results and galaxies in nature. 

This all said, while our models appear to provide a plausi- 
ble match to observed data, the comparisons made thus far are in 
essence postdictions. An even more powerful test of any theoretical 
model would be testable predictions of future surveys. In this sec- 
tion, we sketch out potential observable tests of these models which 

show higher excitation conditions. This lends itself to a direct testable pre- 
diction from these models which we outline in § 17.11 



are imminently possible with the latest generation of bolometer and 
heterodyne receiver technology. 

7.2 CO Linewidths of z ~ 2 SMGs and Quasars 

Our model for the increased line widths from SMGs relies on a tem- 
porary increase in velocity dispersions owing to multiple galaxies 
in the simulation box/telescope beam. Two features of this model 
are evident. First, the line widths are relatively low, and representa- 
tive of the virial velocity of a single progenitor galaxy during the in- 
spiral phase when the observed 850 |J.m flux may be relatively low 
(S l 850~ 1 mJy; Figure[3}. They increase concomitant to the increase 
in submillimeter flux, and are thus broader when the galaxy is in its 
transient SMG phase, though remain broad in the post-SMG phase 
when the galaxies 850 (im flux decreases again. In this picture, one 
might expect a broad range of linewidths from galaxies with lower 
(~1 mJy) 850 \im fluxes: small line widths corresponding to inspi- 
ralling galaxies, and large line widths corresponding to post-SMG 
phase galaxies. Similarly, lower mass mergers (e.g. model SMG 13) 
at final coalescence may contribute to the dispersion of line widths 
on the low flux end. 

Second, as discussed bv lNaravanan et al] |2009), the 850 irm 
lightcurve from merging galaxies has a similar shape for a mass 
sequence of mergers, though scales in normalization. This means 
that the submillimeter flux curve shown in the top panel of Fig- 
ure [3] simply scales upward with increasing merger mass (indeed 
this is evident with a direct compariso n of the lightcurve pre sented 
in Figure [3] with that of Figure 1 in iNaravanan e t al. 2009). This 
means that in more massive mergers, the inspiral phase corresponds 
to 5 mJy sources, and the peak SMG phase to rarer, more luminous 
(S l 850^ 10-15 mJy) SMGs. Consequently, there may be a spread 
in CO line widths from even 5 mJy SMGs. The most luminous 
sources, however, only occur at the final coalescence burst of ex- 
tremely rare ~10 13 Mq mergers. As such, S$5o~ 10-20 mJy galax- 
ies will have extremely broad (FWHM > 1000 km s _1 ) CO line 
widfh^f] with a relatively small dispersion. 

The aforementioned effects will have the following generic 
consequence for CO linewidths from high redshift mergers: the 
mean CO FWHM will increase with increasing submillimeter flux, 
and the dispersion will decrease. We show this explicitly in Fig- 
ure [T2] by plotting the sightline-averaged model CO FWHM ver- 
sus 850 am flux for model galaxies SMG1, SMG 10 and SMGlfl 
in Table |2~T1 for a flux limit of Ss5o>l mJy, and in the middle 
panel, the CO FWHM- 1.1 mm relation. This serves as a predic- 
tion for both future sensitive 85 (im surveys, as well as 1.1 mm 
(e.g. AzTEC; lwilson et alj200^) counts. 

In the right panel of Figure [12] we attempt to compare our 
model pred iction in the left panel w i th dat a culled from the recent 
review by ISolomon & Vanden Bou 3 J2005h . An important point is 
that optically selected quasars cannot be included in this sort of 
comparison as they may preferentially have their molecular disks 
in a face-on configuration, thus skewing avera ge line widths (e.g. 
ICarilli & Wan3l2006l : INaravanan etafll2008ch . This does exclude 
a number of sources which may be si multaneously quasars , and 
submillimeter bright (e.g. Figure[3]and ICoppin et alj|2008b1) . The 



8 Recall that the CO line width scales with circular velocity, and conse- 
quently, with galaxy mass. 

9 These models were chosen to span a broad mass range. Note that model 
SMG 13 is too l ow in mass to form a detectable S$rM>5 rnjy SMG; 
INaravanan et al l 120091) . 
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Figure 12. Predicted CO (J=3-2) F WHM-5 8 50 relation (left p a nel), p redicted CO (J=3-2) FWHM-5i.i mm relation (middle panel), and observed FWHM- 
5 8 50 relation culled from data from lSolomon & Vanden Bouj J2005t> - The models are for models SMG1, SMG10 and SMG13. SMG1 and SMG10 were 
chosen as usual to bracket the range of halo masses employed here, and SMG13 added to probe lower fluxes (though note it is too low in mass to ever form a 
S'850> 5 mly SMG). The brightest SMGs at 850 |om and 1.1 mm are predicted to have a narrow dispersion of CO FWHMs, and typically broad lines whereas 
lower luminosity objects may have a larger range of line widths. The mean CO line width will increase with increasing (sub)mm flux. See text for reasoning. 
From the observational dataset, optical quasars have been eliminated. The comparison to the observations is gener ally inconclusive. This may ow e to many 
uncertain gravitational lens magnifications in the observed SMGs (known lensed sources are marked as red squares; ISolomon & Vanden Boum2005I) . Because 
these models probe both the full dynamic range of simulated 850 (am fluxes and CO FWHMs, adding more models does not increase the dispersion in the 
models plots. 



observed data in the right panel of Figure [T_2] are generally in- 
conclusive in comparison to the predictions in the left panel. The 
effects of uncertain gravitational lens magnifications muddy in- 
terpretation (sources that are known to be lensed are marked as 
red squares; ISolomon & Vanden Boutl 2005). Better statistics with 
SCUBA2 and AzTEC will provide a direct and imminent test for 
this aspect of these models. 

7.3 Predicted CO Excitation for High-Resolution 
Observations of SMGs 

With the advent of broadband receivers on a multitude of tele- 
scopes, CO SEDs of high redshift objects will soon become avail- 
able and serve as a test for the excitation predictions outlined in 

We first note that Figure [TJJ itself is a testable prediction. As 
few multi-line surveys are currently published, Figure [TTI serves as 
a direct comparative for future constraints on the CO SED from 
SMGs. A crucial component to this will be the direct detection of 
CO (J=l-0) emission from SMGs. In the absence of this, relative 
intensities to this transition will be difficult to interpret, and direct 
conversion of emission from higher lying transitions to H2 masses 
unclear. This is currently feasible at a variety of telescope facilities, 
though interferometers (e.g. the EVLA) may be pref erable as base- 
lines may be problematic on single dish telescopes (Hainli ne et al.1 
l200fj) . 

Second, we may make predictions for imminent higher spa- 
tial resolution observations (e.g. ALMA). The CO excitation con- 



ditions presented in Figure [TT] for our model SMGs were aver- 
aged over the entire 8 kpc box, thus including emission from sub- 
thermally excited gas. While the gas associated with the nuclear 
starburst is indeed warm and thermalized, including diffuse, sub- 
thermal gas had the effect of lowering the mean excitation con- 
dition. This is analogous to observations of z ~ 2 SMGs which 
include both dense, nuclear gas, as well as diffuse gas in the tele- 
scope beam. Consequently, the modeled excitation conditions in 
Figure QT] showed excellent agreement with multi-line CO mea- 
surements from SMGs. 



In principle, one can imagine that higher spatial resolution 
observations of SMGs which focus on the dense, nuclear regions 
would show higher excitation CO SEDs. Our approach allows us to 
quantify this effect, as well as make predictions for the next gen- 
eration of high spatial resolution interferometers (e.g. ALMA). In 
Figure[J_3] we plot the sightline-averaged CO SED for SMG models 
SMG1 and SMG 10 as a function of decreasing physical beam size. 
As before, we only consider SMGs with a fiducial selection criteria 
of Ssso^ 5 mJy. As the observations probe narrower and narrower 
columns, the observed gas is on average more highly excited, and 
the turnover point of the CO SED moves to increasingly high ro- 
tational number. The lack of complete thermalization in the lower 
observed transitions (e.g. CO J=2-l) owes to the fact that diffuse 
gas is still folded into the observation, even for rather high spatial 
resolution observations. 
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CO Upper Rotational Quantum Number (J upper ) 

Figure 13. Predicted CO line SEDs for models SMG1 and SMG10 dur- 
ing their SMG phases at increasing spatial resolution. The black solid 
line shows the predicted CO line SED for thermalized populations. As the 
warmer, denser gas toward the nucleus is probed, the peak and shape of the 
excitation ladder shift toward higher levels as more of the gas becomes ther- 
malized. High resolution observations by ALMA of the nucleus of SMGs 
will typically show higher excitation CO line SEDs than lower resolution 
observations. 



8 SUMMARY 

We have combined hydrodynamic simulations of Submillimeter 
Galaxy formation and evolution with 3D non-LTE molecular line 
radiative transfer calculations to provide a model for the CO emis- 
sion properties from SMGs. Our model has shown a number of 
successes in matching the observed spatial extent of CO emission, 
CO line widths and excitation conditions. We utilized these models 
to understand the origin of these emission properties: 

• In our model, SMGs originate in major mergers 
dNaravanan et al.l l2009h . Strong gaseous inflows drive highly 
concentrated molecular gas complexes such that the observed 
characteristic CO radius is of order ~1.5 kpc. The large radius 
tail of the distribution arises from pre-coalescence galaxies in 
extremely massive (~1O 13 M ) halos which are SMGs even 
during the inspiral phase (e.g. Figure 1 of iNaravanan etaiT2 009). 

• The large CO line widths from SMGs owe to the fact that they 
are typically being observed during a transient phase where the gas 
is highly non-virialized and multiple galaxies are in the simula- 
tion box/telescope beam. During interactions, the CO FWHM from 
mergers is roughly 2.35 x the circular velocity of a progenitor (for 
a Gaussian line; Equation[4]l. Two merging V c « 320 km s' 1 disks 
naturally produce average (Ssbo~5-7 mjy) SMGs with extremely 
broad line widths of order ~600-800 km s _1 . 

We have additionally been able to provide interpretation re- 
garding the usage of CO as a diagnostic of physical conditions: 

• The usage of CO line widths from SMGs as a dynamical mass 
estimator may overestimate the enclosed mass. Typical overesti- 
mates are of order Md yn / M ac tuai ~ 1.5-2. 

• The CO excitation in SMGs is high, with the rotational lad- 
der turning over at the ~J=5 or 6 level. The level populations at 
Jupper > 2 are typically subthermal, and assumptions regarding 
brightness temperature ratios of unity with the ground state will 
lead to underestimates of the inferred H2 gas mass. The CO (J=3- 



2) line from SMGs is typically a factor of ~3 below the intensity 
expected from thermalized level populations. 

Finally, we have made predictions for this model which are 
imminently testable with the newest generation of bolometer arrays 
and wide-band sensitive CO receivers: 

• The CO line widths from galaxies which will become SMGs 
are predicted to be their broadest when the submillimeter flux is 
the highest. Consequently, the brightest SMGs (5*850 > 15 mJy) 
are predicted to only have quite broad (FWHM> 1000 km s _1 ) 
line widths with a small dispersion in FWHMs observed. Lower 
flux galaxies (5 , 85o~l) can be either normal disks (where the line 
widths are predicted to be relatively narrow), or post-SMG phase 
mergers (where the line widths will be relatively broad). As such, 
lower flux galaxies are predicted to have a broad dispersion in ob- 
served CO FWHM. Care must be taken to both remove optically 
selected quasars (which may be biased to have face-on molecular 
disks) and lens-magnified sources from the sample. 

• The intensity from the CO J=5-4 or 6-5 line, where the CO 
SED turns over, is predicted to be a factor of ~5-7 times that of 
the ground state. While detections at these higher lying lines are 
becoming routine, future detections of CO (J=l-0) emission from 
SMGs will be required to test the predicted rotational ladders. 

• The peak and shape of the CO SED will shift toward higher 
lying transitions as smaller spatial extents are investigated with fu- 
ture high resolution arrays (e.g. ALMA). 
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